Impacts of longitudinal water curtain cooling system on transcriptome-related immunity in ducks.

BACKGROUND: The closed poultry houses integrated with a longitudinal water curtain cooling system (LWCCS) are widely used in modern poultry production. This study showed the variations in environmental conditions in closed houses integrated with a longitudinal water curtain cooling system. We evaluated the influence of different environmental conditions on duck growth performance and the transcriptome changes of immune organs, including the bursa of Fabricius and the spleen. RESULT: This study investigated the slaughter indicators and immune organ transcriptomes of 52-day-old Cherry Valley ducks by analyzing the LWCC at different locations (water curtain end, middle position, and fan cooling end). The results showed that the cooling effect of the LWCCS was more evident from 10:00 a.m. -14:00. And from the water curtain end to the fan cooling end, the hourly average temperature differently decreased by 0.310℃, 0.450℃, 0.480℃, 0.520℃, and 0.410℃, respectively (P < 0.05). The daily and hourly average relative humidity decreased from the water curtain end to the fan cooling end, dropping by 7.500% and 8.200%, respectively (P < 0.01). We also observed differences in production performance, such as dressing weight, half-eviscerated weight, skin fat rate, and percentage of abdominal fat (P < 0.01), which may have been caused by environmental conditions. RNA-sequencing (RNA-seq) revealed 211 and 279 differentially expressed genes (DEGs) in the ducks' bursa of Fabricius and spleen compared between the water curtain end and fan cooling end, respectively. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the two organs showed the DEGs were mainly enriched in cytokine-cytokine receptor interaction, integral component of membrane, Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) signaling pathway, etc. Our results implied that full-closed poultry houses integrated with LWCCS could potentially alter micro-environments (water curtain vs. fan cooling), resulting in ducks experiencing various stressful situations that eventually affect their immunity and production performance. CONCLUSION: In this study, our results indicated that uneven distributions of longitudinal environmental factors caused by LWCCS would affect the dressed weight, breast muscle weight, skin fat rate, and other product performance. Moreover, the expression of immune-related genes in the spleen and bursa of ducks could be affected by the LWCCS. This provides a new reference to optimize the use of LWCCS in conjunction with close duck houses in practical production.

Genome-wide association study for bone quality of ducks during the laying period.

The cage-rearing model of the modern poultry industry makes the bones of birds, especially egg-laying birds, more vulnerable to fracture, which poses serious damage to the health of birds. Research confirms that genetic material plays an important role in regulating bone growth, development, and remodeling. However, the genetic architecture underlying bone traits is not well understood. The objectives of this study are to identify valuable genes and genetic markers through a genome-wide association study (GWAS) for breeding to improve the duck bone quality. First, we quantified the tibia and femur quality traits of 260 laying ducks. Based on GWAS, a total of 75 SNP loci significantly associated with bone quality traits were identified, and 67 potential candidate genes were annotated. According to gene function analysis, genes P4HA2, WNT3A, and BST1 et al may influence bone quality by regulating bone cell activity, calcium and phosphate metabolism, or bone collagen maturation and cross-linking. Meanwhile, combined with the transcriptome results, we found that HOXB cluster genes are also important in bone growth and development. Therefore, our findings were helpful in further understanding the genetic architecture of the duck bone quality and provided a worthy theoretical basis and technological support to improve duck bone quality by breeding.

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes.

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.

Identification of candidate genes affecting the tibia quality in Nonghua duck.

The skeleton is a vital organ providing structural support in poultry. Weakness in bone structure can lead to deformities, osteoporosis, cage fatigue, and fractures, resulting in economic losses. Research has substantiated that genetic factors play a significant role in influencing bone quality. The discovery of genetic markers associated with bone quality holds paramount importance for enhancing genetic traits related to the skeletal system in poultry. This study analyzed nine phenotypic indicators of tibia quality in 120-day-old ducks. The phenotypic correlation revealed a high correlation among diameter, Perimeter, and weight (0.69-0.78), and a strong correlation was observed between toughness and breaking strength (0.62). Then, we conducted a genome-wide association analysis of the phenotypic indicators to elucidate the genetic basis of tibial quality in Nonghua ducks. Among the 11 candidate genes that were annotated, TAPT1, BST1, and STIM2 were related to the diameter indicator, ZNF652, IGF2BP1, CASK, and GREB1L were associated with the weight and toughness indicators. RFX8, GLP1R, and DNAAF5 were identified for ash, calcium, and phosphorus content, respectively. Finally, KEGG and GO analysis for annotated genes were performed. STIM2 and BST1 were enriched into the Calcium signalling pathway and Niacin and nicotinamide metabolic pathway, which may be key candidate genes affecting bone quality phenotypes. Gene expression analysis of the candidate genes, such as STIM2, BST1, TAPT1, and CASK showed higher expression levels in bones compared to other tissues. The obtained results can contribute to new insights into tibial quality and provide new genetic biomarkers that can be employed in duck breeding.

Genome sequencing of drake semen micobiome with correlation with their compositions, sources and potential mechanisms affecting semen quality.

Artificial insemination (AI) technology has greatly promoted the development of the chicken industry. Recently, AI technology has also begun to be used in the duck industry, but there are some problems. Numerous researchers have shown that microbes colonizing in semen can degrade semen quality, and AI can increase the harmful microbial load in hen's reproductive tract. Different from the degraded external genitalia of roosters, drakes have well-developed external genitalia, which may cause drake semen to be more susceptible to microbial contamination. However, information on the compositions, sources, and effects of semen microbes on semen quality remains unknown in drakes. In the current study, high-throughput sequencing technology was used to detect microbial communities in drake semen, environmental swabs, cloacal swabs, and the spermaduct after quantifying the semen quality of drakes to investigate the effects of microbes in the environment, cloaca, and spermaduct on semen microbiota and the relationships between semen microbes and semen quality. Taxonomic analysis showed that the microbes in the semen, environment, cloaca, and spermaduct samples were all classified into 4 phyla and 25 genera. Firmicutes and Proteobacteria were the dominant phyla. Phyllobacterium only existed in the environment, while Marinococcus did not exist in the cloaca. Of the 24 genera present in semen: Brachybacterium, Brochothrix, Chryseobacterium, Kocuria, Marinococcus, Micrococcus, Rothia, Salinicoccus, and Staphylococcus originated from the environment; Achromobacter, Aerococcus, Corynebacterium, Desemzia, Enterococcus, Jeotgalicoccus, Pseudomonas, Psychrobacter, and Turicibacter originated from the cloaca; and Agrobacterium, Carnobacterium, Chelativorans, Devosia, Halomonas, and Oceanicaulis originated from the spermaduct. In addition, K-means clustering analysis showed that semen samples could be divided into 2 clusters based on microbial compositions, and compared with cluster 1, the counts of Chelativorans (P < 0.05), Devosia (P < 0.01), Halomonas (P < 0.05), and Oceanicaulis (P < 0.05) were higher in cluster 2, while the sperm viability (P < 0.05), total sperm number (P < 0.01), and semen quality factor (SQF) (P < 0.01) were lower in cluster 2. Furthermore, functional prediction analysis of microbes showed that the activities of starch and sucrose metabolism, phosphotransferase system, ABC transporters, microbial metabolism in diverse environments, and quorum sensing pathways between cluster 1 and cluster 2 were significantly different (P < 0.05). Overall, environmental/cloacal microbes resulted in semen contamination, and microbes from the Chelativorans, Devosia, Halomonas, and Oceanicaulis genera may have negative effects on semen quality in drakes by affecting the activities of starch and sucrose metabolism, phosphotransferase system, ABC transporters, and quorum sensing pathways that are associated with carbohydrate metabolism. These data will provide a basis for developing strategies to prevent microbial contamination of drake semen.

Multi-omics reveals goose fatty liver formation from metabolic reprogramming.

To comprehensively provide insight into goose fatty liver formation, we performed an integrative analysis of the liver transcriptome, lipidome, and amino acid metabolome, as well as peripheral adipose tissue transcriptome analysis using samples collected from the overfed geese and normally fed geese. Transcriptome analysis showed that liver metabolism pathways were mainly enriched in glucolipid metabolism, amino acid metabolism, inflammation response, and cell cycle; peripheral adipose tissue and the liver cooperatively regulated liver lipid accumulation during overfeeding. Liver lipidome patterns obviously changed after overfeeding, and 157 different lipids were yielded. In the liver amino acid metabolome, the level of Lys increased after overfeeding. In summary, this is the first study describing goose fatty liver formation from an integrative analysis of transcriptome, lipidome, and amino acid metabolome, which will provide a whole new dimension to understanding the mechanism of goose fatty liver formation.

Cascade Co8FeS8@Co1-xS nano-enzymes trigger efficiently apoptosis-ferroptosis combination tumor therapy.

This study involved the preparation of Metal Organic Frameworks (MOF)-derived Co8FeS8@Co1-xS nanoenzymes with strong interfacial interactions. The nanoenzymes presented the peroxidase (POD)-like activity and the oxidation activity of reduced glutathione (GSH). Accordingly, the dual activities of Co8FeS8@Co1-xS provided a self-cascading platform for producing significant amounts of hydroxyl radical (•OH) and depleting reduced glutathione, thereby inducing tumor cell apoptosis and ferroptosis. More importantly, the Co8FeS8@Co1-xS inhibited the anti-apoptosis protein B-cell lymphoma-2 (Bcl-2) and activated caspase family proteins, which caused tumor cell apoptosis. Simultaneously, Co8FeS8@Co1-xS affected the iron metabolism-related genes such as Heme oxygenase-1 (Hmox-1), amplifying the Fenton response and promoting apoptosis and ferroptosis. Therefore, the nanoenzyme synergistically killed anti-apoptotic tumor cells carrying Kirsten rat sarcoma viral oncogene homolog (KRAS) mutations. Furthermore, Co8FeS8@Co1-xS demonstrated good biocompatibility, which paved the way for constructing a synergistic catalytic nanoplatform for an efficient tumor treatment.

Molecular mechanisms of hypothalamic-pituitary-ovarian/thyroid axis regulating age at first egg in geese.

Age at first egg (AFE) has consistently garnered interest as a crucial reproductive indicator within poultry production. Previous studies have elucidated the involvement of the hypothalamic-pituitary-ovarian (HPO) and hypothalamic-pituitary-thyroid (HPT) axes in regulating poultry sexual maturity. Concurrently, there was evidence suggesting a potential co-regulatory relationship between these 2 axes. However, as of now, no comprehensive exploration of the key pathways and genes responsible for the crosstalk between the HPO and HPT axes in the regulation of AFE has been reported. In this study, we conducted a comparative analysis of morphological differences and performed transcriptomic analysis on the hypothalamus, pituitary, thyroid, and ovarian stroma between normal laying group (NG) and abnormal laying group (AG). Morphological results showed that the thyroid index difference (D-) value (thyroid index D-value=right thyroid index-left thyroid index) was significantly (P < 0.05) lower in the NG than in the AG, while the ovarian index was significantly (P < 0.01) higher in the NG than in the AG. Furthermore, between NG and AG, we identified 99, 415, 167, and 1182 differentially expressed genes (DEGs) in the hypothalamus, pituitary, thyroid, and ovarian stroma, respectively. Gene ontology (GO) analysis highlighted that DEGs from 4 tissues were predominantly enriched in the "biological processes" category. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that 16, 14, 3, and 26 KEGG pathways were significantly enriched (P < 0.05) in the hypothalamus, pituitary, thyroid, and ovarian stroma. The MAPK signaling pathway emerged as the sole enriched pathway across all 4 tissues. Employing an integrated analysis of the protein-protein interaction (PPI) network and correlation analysis, we found GREB1 emerged as a pivotal component within the HPO axis to regulate estrogen-related signaling in the HPT axis, meanwhile, the HPT axis influenced ovarian development by regulating thyroid hormone-related signaling mainly through OPN5. Then, 10 potential candidate genes were identified, namely IGF1, JUN, ERBB4, KDR, PGF, FGFR1, GREB1, OPN5, DIO3, and THRB. These findings establish a foundation for elucidating the physiological and genetic mechanisms by which the HPO and HPT axes co-regulate goose AFE.

Transcriptome-based comparison reveals key genes regulating allometry growth of forelimb and hindlimb bone in duck embryos.

Allometric growth of the forelimb and hindlimb is a widespread phenomenon observed in vertebrates. As a typical precocial bird, ducks exhibit more advanced development of their hindlimbs compared to their forelimbs, enabling them to walk shortly after hatching. This phenomenon is closely associated with the development of long bones in the embryonic stage. However, the molecular mechanism governing the allometric growth of duck forelimb and hindlimb bones is remains elusive. In this study, we employed phenotypic, histological, and gene expression analyses to investigate developmental differences between the humerus (forelimb bone) and tibia/femur (hindlimb bones) in duck embryos. Our results revealed a gradual increase in weight and length disparity between the tibia and humerus from E12 to E28 (embryo age). At E12, endochondral ossification was observed solely in the tibia but not in the humerus. The number of differentially expressed genes (DEGs) gradually increased at H12 vs. T12, H20 vs. T20, and H28 vs. T28 stages consistent with phenotypic variations. A total of 38 DEGs were found across all 3 stages. Protein-protein interaction network analysis demonstrated strong interactions among members of HOXD gene family (HOXD3/8/9/10/11/12), HOXB gene family (HOXB8/9), TBX gene family (TBX4/5/20), HOXA11, SHOX2, and MEIS2. Gene expression profiling indicated higher expression levels for all HOXD genes in the humerus compared to tibia while opposite trends were observed for HOXA/HOXB genes with low or no expression detected in the humerus. These findings suggest distinct roles played by different clusters within HOX gene family during skeletal development regulation of duck embryo's forelimbs versus hind limbs. Notably, TBX4 exhibited high expression levels specifically in tibia whereas TBX5 showed similar patterns exclusively within humerus as seen previously across other species' studies. In summary, this study identified key regulatory genes involved in allometric growth of duck forelimb and hindlimb bones during embryonic development. Skeletal development is a complex physiological process, and further research is needed to elucidate the regulatory role of candidate genes in endochondral ossification.

Multi-Enzyme Activity of MIL-101 (Fe)-Derived Cascade Nano-Enzymes for Antitumor and Antimicrobial Therapy.

The clinical application of oncology therapy is hampered by high glutathione concentrations, hypoxia, and inefficient activation of cell death mechanisms in cancer cells. In this study, Fe and Mo bimetallic sulfide nanomaterial (FeS2 @MoS2 ) based on metal-organic framework structure is rationally prepared with peroxidase (POD)-, catalase (CAT)-, superoxide dismutase (SOD)-like activities and glutathione depletion ability, which can confer versatility for treating tumors and mending wounds. In the lesion area, FeS2 @MoS2 with SOD-like activity can facilitate the transformation of superoxide anions (O2 - ) to hydrogen peroxide (H2 O2 ), and then the resulting H2 O2 serves as a substrate for the Fenton reaction with FMS to produce highly toxic hydroxyl radicals (∙OH). Simultaneously, FeS2 @MoS2 has an ability to deplete glutathione (GSH) and catalyze the decomposition of nicotinamide adenine dinucleotide phosphate (NADPH) to curb the regeneration of GSH from the source. Thus it can realize effective tumor elimination through synergistic apoptosis-ferroptosis strategy. Based on the alteration of the H2 O2 system, free radical production, glutathione depletion and the alleviation of hypoxia in the tumor microenvironment, FeS2 @MoS2 NPS can not only significantly inhibit tumors in vivo and in vitro, but also inhibit multidrug-resistant bacteria and hasten wound healing. It may open the door to the development of cascade nanoplatforms for effective tumor treatment and overcoming wound infection.

Comparative transcriptome analysis reveals mechanisms of restriction feeding on lipid metabolism in ducks.

Presently, excessive fat deposition is the main reason to limit the development of duck industry. In the production, the methods of restricted feeding (RF) were widely used to reduce the lipid deposition of ducks. The liver (L), abdominal adipose (AA), and subcutaneous adipose (SA) were the main tissues of lipid metabolism and deposition of ducks. However, the mechanisms of lipid metabolism and deposition of ducks under RF have not been fully clarified. In this study, in order to better understand the mechanisms of lipid metabolism and deposition in ducks under RF, a total of 120 male Nonghua ducks were randomly divided into a free feeding group (FF, n = 60) and RF group (RF, n = 60), then comparative transcriptomic analysis of L, AA, and SA between FF (n = 3) and RF (n = 3) ducks was performed at 56 d of age. Phenotypically, L, AA, and SA index of FF group was higher than that in RF group. There were 279, 390, and 557 differentially expressed genes (DEGs) in L, AA, and SA. Functional enrichment analysis revealed that ECM-receptor interaction and metabolic pathways were significantly enriched in L, AA, and SA. Lipid metabolism-related pathways including fatty acid metabolism, unsaturated fatty acid synthesis, and steroidogenesis were significantly enriched in AA and SA. Moreover, through integrated analysis weighted gene coexpression network (WGCNA) and protein-protein interaction network, 10 potential candidate genes involved in the ECM-receptor interaction and lipid metabolism pathways were identified, including 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2), aldolase B (ALDOB), formimidoyltransferase cyclodeaminase(FTCD), phosphoenolpyruvate carboxykinase 1 (PCK1), tyrosine aminotransferase (TAT), stearoyl-CoA desaturase (SCD), squalene epoxidase (SQLE), phosphodiesterase 4B (PDE4B), choline kinase A (CHKA), and elongation of very-long-chain fatty acids-like 2 (ELOVL2), which could play a key role in lipid metabolism and deposition of ducks under RF. Our study reveals that the liver might regulate the lipid metabolism of abdominal adipose and subcutaneous adipose through ECM-receptor interaction and metabolic pathways (fatty acid metabolism, unsaturated fatty acid synthesis, and steroid synthesis), thus to reduce the lipid deposition of ducks under RF. These results provide novel insights into the avian lipid metabolism and will help better understand the underlying molecular mechanisms.

Comprehensive transcriptomic and metabolomic analysis of the effect of feed restriction on duck sternal development.

Skeletal characteristics are important to the growth and development of poultry. In feeding management, constant free feeding (FF) of poultry may lead to imbalance between bone development and weight gain. Feed restriction (FR), to a certain extent, is one way to solve this problem. However, the effect of feed restriction on poultry bone development needs further elucidation at the molecular level. Therefore, in the present study, we investigated the effects of different levels of feed restriction (60% FR, 70% FR, 80% FR, and FF) on the sternum development of ducks at 7 and 8 wk old. In the seventh wk, with increasing feed restriction, the values of traits including body weight, breast muscle weight, sternal weight, keel length, and calcified keel length decreased. However, in the eighth wk, the sternum weight and keel length of ducks treated with 60% FR were unexpectedly higher than those of FF individuals, indicative of catch-up growth. Then, we conducted RNA-seq and metabolomic analysis on sterna from 7- and 8-wk-old FF and 60% FR ducks. The results identified multiple differentially expressed genes (DEGs) associated with sternum development that were influenced by feed restriction. Among them, we found that the mRNA expression levels of the chondroitin sulfate synthase 3 (CHSY3) and annexin A2 (ANXA2) which are involved in glycosaminoglycan biosynthesis and bone mineralization, had smaller changes over time under FR treatment than under FF treatment, implying that the FR treatment to a certain extent prevented the premature calcification and prolonged the development time of duck sternum. In addition, the metabolomic and integrative analyses revealed that several antiaging-related metabolites and genes were associated with sternal catch-up growth. Pyrimidine metabolism was identified as the most significant pathway in which most differential metabolites (DMs) between FF and 60% FR were enriched. The results from integrative analysis revealed that the content and expression of 4-aminobutyric acid (GABA) and its related genes showed relatively higher activity in the 60% FR group than in the FF group. The present study identifies multiple biomarkers associated with duck sternum development that are influenced by feed restriction and suggests the potential mechanism of feed restriction-associated duck sternal catch-up growth.

Estrogen Receptor Gene 1 (ESR1) Mediates Lipid Metabolism in Goose Hierarchical Granulosa Cells Rather than in Pre-Hierarchical Granulosa Cells.

(1) Background: The role of estrogen receptor gene 1 (ESR1) in female reproduction and lipid metabolism has been extensively investigated. However, its contribution to lipid metabolism during the development of poultry follicles remains unclear. (2) Methods: This study aimed to explore the function of ESR1 via overexpressing (ESR1ov) and interfering (ESR1si) with its expression in pre-hierarchical granulosa cells (phGCs) and hierarchical granulosa cells (poGCs). (3) Results: We successfully cloned and obtained an 1866 bp segment of the full-length CDS region of the Sichuan white goose ESR1 gene. In phGCs of the ESR1ov and ESR1si groups, there were no significant changes compared to the control group. However, in poGCs, the ESR1ov group exhibited decreased lipid deposition, triglycerides, and cholesterol compared to the control group, while the ESR1si group showed increased lipid deposition, triglycerides, and cholesterol. The expression of APOB and PPARα was significantly reduced in the ESR1ov group compared to the ESR1ov-NC group. Moreover, significant changes in the expression of ACCα, DGAT1, SCD, CPT1, and ATGL were observed between the ESR1si and ESR1si-NC group. (4) Conclusions: These findings shed light on the function and molecular mechanism of ESR1 in lipid metabolism in goose poGCs, providing a better understanding of the physiological process of goose follicular development.

Exploring right ovary degeneration in duck and goose embryos by histology and transcriptome dynamics analysis.

BACKGROUND: The development of asymmetric chick gonads involves separate developmental programs in the left and right gonads. In contrast to the left ovary developing into a fully functional reproductive organ, the right ovary undergoes gradual degeneration. However, the molecular mechanisms underlying the the degeneration of the right ovary remain incompletely understood. In the present study, we investigated the histomorphological and transcriptomic changes in the right ovary of ducks and geese during the the embryonic stage up to post-hatching day 1. RESULT: Hematoxylin-eosin stainings revealed that the right ovary developed until embryonic day 20 in ducks (DE20) or embryonic day 22 in geese (GE22), after which it started to regress. Further RNA-seq analyses revealed that both the differentially expressed genes (DEGs) in ducks and geese right ovary developmental stage were significantly enriched in cell adhesion-related pathway (ECM-receptor interaction, Focal adhesion pathway) and Cellular senescence pathway. Then during the degeneration stage, the DEGs were primarily enriched in pathways associated with inflammation, including Herpes simplex virus 1 infection, Influenza A, and Toll-like receptor signaling pathway. Moreover, duck-specific DEGs showed enrichment in Steroid hormone biosynthesis, Base excision repair, and the Wnt signaling pathway, while geese-specifically DEGs were found to be enriched in apoptosis and inflammation-related pathways, such as Ferroptosis, Necroptosis, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway. These findings suggest that the degeneration process of the right ovary in ducks occurs at a slower pace compared to that in geese. Additionally, the observation of the left ovary of the geese varying degeneration rates in the right ovary after hatching indicated that the development of the left ovary may be influenced by the degeneration of the right ovary. CONCLUSION: The data presented in this study provide valuable insights into the dynamic changes in histological structure and transcriptome during the degeneration of the right ovary in ducks and geese. In addition, through the analysis of shared characteristics in the degeneration process of the right ovary in both ducks and geese, we have uncovered the patterns of degradation and elucidated the molecular mechanisms involved in the regression of the right ovary in poultry. Furthermore, we have also made initial discoveries regarding the relationship between the degeneration of the right ovary and the development of the left ovary.

Integration of GWAS and eGWAS to screen candidate genes underlying green head traits in male ducks.

Sexually dimorphic plumage coloration is widespread in birds. The male possesses more brightly colored feathers than the female. Dark green head feathers comprise one of the most typical appearance characteristics of the male Ma duck compared with the female. However, there are noticeable individual differences observed in these characteristics. Herein, genome-wide association studies (GWAS) were employed to investigate the genetic basis of individual differences in male duck green head-related traits. Our results showed that 165 significant SNPs were associated with green head traits. Meanwhile, 71 candidate genes were detected near the significant SNPs, including four genes (CACNA1I, WDR59, GNAO1 and CACNA2D4) related to the individual differences in the green head traits of male ducks. Additionally, the eGWAS identified three SNPs located within two candidate genes (LOC101800026 and SYNPO2) associated with TYRP1 gene expression, and might be important regulators affecting the expression level of TYRP1 in the head skin of male ducks. Our data also suggested that transcription factor MXI1 might regulate the expression of TYRP1, thereby causing differences in the green head traits among male ducks. This study provided primary data for further analysis of the genetic regulation of duck feather color.

Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis.

BACKGROUND: The genetic locus responsible for duck body size has been fully explained before, but the growth trait-related genetic basis is still waiting to be explored. For example, the genetic site related to growth rate, an important economic trait affecting marketing weight and feeding cost, is still unclear. Here, we performed genome wide association study (GWAS) to identify growth rate-associated genes and mutations. RESULT: In the current study, the body weight data of 358 ducks were recorded every 10 days from hatching to 120 days of age. According to the growth curve, we evaluated the relative and absolute growth rates (RGR and AGR) of 5 stages during the early rapid growth period. GWAS results for RGRs identified 31 significant SNPs on autosomes, and these SNPs were annotated by 24 protein-coding genes. Fourteen autosomal SNPs were significantly associated with AGRs. In addition, 4 shared significant SNPs were identified as having an association with both AGR and RGR, which were Chr2: 11483045 C>T, Chr2: 13750217 G>A, Chr2: 42508231 G>A and Chr2: 43644612 C>T. Among them, Chr2: 11483045 C>T, Chr2: 42508231 G>A, and Chr2: 43644612 C>T were annotated by ASAP1, LYN and CABYR, respectively. ASAP1 and LYN have already been proven to play roles in the growth and development of other species. In addition, we genotyped every duck using the most significant SNP (Chr2: 42508231 G>A) and compared the growth rate difference among each genotype population. The results showed that the growth rates of individuals carrying the Chr2: 42508231 A allele were significantly lower than those without this allele. Moreover, the results of the Mendelian randomization (MR) analysis supported the idea that the growth rate and birth weight had a causal effect on the adult body weight, with the growth rate having a greater effect size. CONCLUSION: In this study, 41 SNPs significantly related to growth rate were identified. In addition, we considered that the ASAP1 and LYN genes are essential candidate genes affecting the duck growth rate. The growth rate also showed the potential to be used as a reliable predictor of adult weight, providing a theoretical reference for preselection.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis.

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

Genetic fine-mapping reveals single nucleotide polymorphism mutations in the MC1R regulatory region associated with duck melanism.

Birds are among the most colourful terrestrial vertebrates, with various plumage colours and patterns. We conducted a genome-wide association study (GWAS) on an intercross F2 population of Pekin ducks and mallards (n = 722) and identified a 1.57-Mb genetic region (Chr11: 20,176,480-21,750,101 bp) related to duck melanism. Fine mapping by linkage disequilibrium (LD) and FST analysis narrowed the final candidate region to a region of 22,500 bp (Chr11: 20,677,500-20,700,000 bp) including three coding genes, TCF25, MC1R and TUBB3. Combined with transcriptome and qRT-PCR analysis, MC1R was identified as the unique genetic locus responsible for black plumage in ducks, and it was significantly more highly expressed in the feather bulbs of black ducks. We also identified 52G > A (Chr11: 20,696,354G > A) and 376G > A (Chr11: 20,696,678G > A) mutations in the MC1R coding region that have been widely studied in ducks. In addition, structural variations (SVs) were screened by nanopore sequencing, and no significant SV was found to be associated with the duck black plumage trait. However, we identified four novel single nucleotide polymorphisms in the MC1R regulator region (Chr11: 20,678,412G > A, Chr11: 20,679,236G > A, Chr11: 20,692,496 A > G and Chr11: 20,692,791 A > G) that had a strong association with the black plumage phenotype of ducks and combined with potential changes in transcription binding affinities. The luciferase reporter gene assay demonstrated that Chr11: 20,678,412G > A and Chr11: 20,679,236G > A led to significant promoter activity changes. Our research emphasizes the importance of MC1R regulatory region mutation in determining the duck black plumage phenotype, and these results expand our understanding of the genetic mechanism underlying duck plumage colour.

Ferritin heavy chain participated in ameliorating 3-nitropropionic acid-induced oxidative stress and apoptosis of goose follicular granulosa cells.

Oxidative stress is the major culprits responsible for ovarian dysfunction by damaging granulosa cells (GCs). Ferritin heavy chain (FHC) may participate in the regulation of ovarian function by mediating GCs apoptosis. However, the specific regulatory function of FHC in follicular GCs remains unclear. Here, 3-nitropropionic acid (3-NPA) was utilized to establish an oxidative stress model of follicular GCs of Sichuan white geese. To explore the regulatory effects of FHC on oxidative stress and apoptosis of primary GCs in geese by interfering or overexpressing FHC gene. After transfection of siRNA-FHC to GCs for 60 h, the expressions of FHC gene and protein decreased significantly (P < 0.05). After FHC overexpression for 72 h, the expressions of FHC mRNA and protein upregulated considerably (P < 0.05). The activity of GCs was impaired after interfering with FHC and 3-NPA coincubated (P < 0.05). When overexpression of FHC combined with 3-NPA treatment, the activity of GCs was remarkably enhanced (P < 0.05). After interference FHC and 3-NPA treatment, NF-κB and NRF2 gene expression decreased (P < 0.05), the intracellular reactive oxygen species (ROS) level increased greatly (P < 0.05), BCL-2 expression reduced, BAX/BCL-2 ratio intensified (P < 0.05), the mitochondrial membrane potential decreased notably (P < 0.05), and the apoptosis rate of GCs aggravated (P < 0.05). While overexpression of FHC combined with 3-NPA treatment could promote BCL-2 protein expression and reduce BAX/BCL-2 ratio, indicating that FHC regulated the mitochondrial membrane potential and apoptosis of GCs by mediating the expression of BCL-2. Taken together, our research manifested that FHC alleviated the inhibitory effect of 3-NPA on the activity of GCs. FHC knockdown could suppress the expression of NRF2 and NF-κB genes, reduce BCL-2 expression and augment BAX/BCL-2 ratio, contributing to the accumulation of ROS and jeopardizing mitochondrial membrane potential, as well as exacerbating GCs apoptosis.

Comparative transcriptome analysis identified crucial genes and pathways affecting sperm motility in the reproductive tract of drakes with different libido.

Libido can affect the semen quality of male, and the sperm motility in semen quality parameters is a reliable index to evaluate the fertility of male. In drakes, the sperm motility is gradually acquired in testis, epididymis, and spermaduct. However, the relationship between libido and sperm motility in drakes has not been reported and the mechanisms of testis, epididymis, and spermaduct regulating the sperm motility of drakes are unclear. Therefore, the purpose of the present study was to compare the semen quality of drakes with libido level 4 (LL4) and libido level 5 (LL5), and tried to identify the mechanisms regulating the sperm motility in drakes by performing RNA-seq in testis, epididymis, and spermaduct. Phenotypically, the sperm motility of drakes (P < 0.01), weight of testis (P < 0.05), and organ index of epididymis (P < 0.05) in the LL5 group were significantly better than those in LL4 group. Moreover, compared with the LL4 group, the ductal square of seminiferous tubule (ST) in testis was significantly bigger in the LL5 group (P < 0.05), and the seminiferous epithelial thickness (P < 0.01) of ST in testis and lumenal diameter (P < 0.05) of ductuli conjugentes/dutus epididymidis in epididymis were significantly longer in the LL5 group. In transcriptional regulation, in addition to KEGG pathways related to metabolism and oxidative phosphorylation, lots of KEGG pathways associated with immunity, proliferation, and signaling were also significantly enriched in testis, epididymis, and spermaduct, respectively. Furthermore, through the integrated analysis of coexpression network and protein-protein interaction network, 3 genes (including COL11A1, COL14A1, and C3AR1) involved in protein digestion and absorption pathway and Staphylococcus aureus infection pathway were identified in testis, 2 genes (including BUB1B and ESPL1) involved in cell cycle pathway were identified in epididymis, and 13 genes (including DNAH1, DNAH3, DNAH7, DNAH10, DNAH12, DNAI1, DNAI2, DNALI1, NTF3, ITGA1, TLR2, RELN, and PAK1) involved in Huntington disease pathway and PI3K-Akt signaling pathway were identified in spermaduct. These genes could play crucial roles in the sperm motility of drakes with different libido, and all data the present study obtained will provide new insights into the molecular mechanisms regulating sperm motility of drakes.